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INTRODUCTION 


The DLP family of transmission line protection systems utilize 
wave-form sampling of the current and voltage inputs and employ 
appropriate algorithme and multiple пісгоргосеввокв to implement the 
protective functions and peripheral functions such as fault 
location. This type of relay is commonly referred to as a "digital" 
or "numerical" relay. 


This document provides additional detailed information about the 
DLE family which may not be found in the various individual 
instruction books. The following material covers DLP revisions A, 
B, and С unless explicitly stated otherwise Іп a particular section. 
Refer to the instruction book of a particular DLP relay to determine 
which functions described below are contained in that relay. 


DICITAL RELAY CONCEPTS - 


The basic concepts of digital relay design and operation are 
discussed in severa] IEEE publications and other technical papers. 
The IEEE tutarial course "Microprocessor Relays and ProtecLiun 
Systems", publicaLion number ВВЕНО269-1-РНЕ, is a good reference, 


Sampling Rate - 


The DLP relays sample current and voltage 16 times per cycle. 
Via a setting, The DLP relays operate at either a nominal 60 Hz or 
50 На system frequency, and "track" frequency excursions 15% ОБЕ 
nominal by adjusting tha sampling interval. For example, the 
sampling interval at 60 Hz is 1.0417 milliseconds, and 1% 18 1.0776 
Milliseconds at 58 Hz. 


Discrete Fourier Transform - 


The DLP relays use a recursive Discrete Fourier Transform (DFT) 
to create phasor quantities from the sampled values of current and 
voltage.  Frior to input into the recursive DFT, the sampled values 
of current are first processed to remove any rc affset that may be 
present, А "software transactor! algorithm is used to remove the dc 
offset. А "transacLor" is a Lransformger incorporating an air gap in 
the magnetic core; current 12 applied to the primary winding and 
voltage is obtained from the secondary winding. Because of the air 
gap, the secondary voltage, v, is related to the primary current, I, 
by the expression: 


V = Та 
whara: 7 = transter impedance of the transactor 


Transactors have been used extensively іп electroamachanical and 
static analog relays to effectively remove dc offset from the 
current. In addition, Lhe Lransfer impedance, 2, establishes the 
base reach magnitude and angle of maximum reach (torque) for 
distance functions. 


The transactor's transfer impedance, 7, is also known as a 
"replica" or "mimic" impedance. if а current із passed through an 
impedance that is a replica of the power system impedance, the 
resulting voltage across the replica impedance will not contain any 
Яс offset that may be present in the current. This is а well known 
technique that is used extensively in relay designs, and the 
transactor 1в one implementation of this technique. 


The general equation for fault current which includes de offset 


is: 
-t/T 
i(t) = I*sin(w*t + а - Ө) - i*sin(« - 8j*e 

where: 1 = реак Value of current 
wom 2ХЯ%Ғ 
а = angle on voltage wave at which fault occurs 
Ө = arctan (w*L/R} = impedance angle of power system 
T= L/R {of power system) 


An impedance, R2 + jw*L2, is considered to Бе a replica impedance 
if: 


L/R — L2/R2 = ог ~ 
а = arctan(w*L/R) = 02 = arctan(w*L2/HR2) 


The voltage across this replica impedance is: 


d 
wit) = i(t)*R2 + — 1(Ъ)#12 
dt 
-t/T 
i(t)*R2 = R2*I*ein(w*t + о = Ө) = HBZ2*I*ein(& - 8)*e 
di R -t/T 
— » L2 = w*l2*I*cos(wt + а ~ Ө) + 12% —— *I*sin(a - O)*e 
dt 1. 
Тһе exponential part of vit) is: 
R -trr 
( LZ* —— - H2 )*1*sin(a - 8)*e 
L 
Since R/L = R2/L2, 
n2 
(12%--- R2) - 0 
L2 


and the dc offset is removed. The voltage is then 


v(t) = R2*I*sin(w*t + œ — Ө) + w*L2*I*cos(w*t + а = Ө) 


Since, 


A*sin(w*L) + B*cus(w*t) = (А: + в)? ста (изъ + arctan(B/A]) 
then 


v (t) [ (R2)? 1 (weL2y? 13 I*sin(w*t ta- B + arctan(w*L2/R2)) 


v(t) = [(R2)'^(w*L2)!]* I*sin(w*t + а) 

The conclusion is that the voltage across a replica impedance does 
not contain any dc offset and it is phase shifted e degrees laading 
compared to the current. Appendix I shows the same results 
graphically using a PC mathematics software package. 


In the DLP relays, all current values are processed by the 
"software transactor" algorithm. This algorithm vses consacutive 
samples of current at "t" and "t + 4" where 4 is the Sampling 
interval. The desired quantity is 15 = (R*i + Ledi/dt} where, 
ideally, L/R in the DLP has the same value as the L/R сЕ the 
transmission line. The DLP calculates L/R from the values of two 
settings POSANG and ZERANG. Consequently, it is lmportant that both 
POSANG and ZERANG be set as слова to their respective actual line 
angles as possible. 


(R*i + L*di/dt) is approximated at time "t + 6/2" which is the 
mid-point between two samples. The approximation is: 


iftté) + ift) i(tt6) - i(t) 
iZ = В сее + L —— 
2 6 


All currents are converted to 12 samples and properly scaled prior 
to being used by the recursive DFT algorithm. 


The "Iull-cycle" ОРТ calculations implemented within the DLP are 
shown below: 


N 
2 ӨТЕ 
REAL = — | Sk sin —— 
М М 
k-1 
М 
2 >. РФК 
IMAG = --- Są COR —— 
H H 
k=1 


where: N = 16 samples 


Sy = magnitude of kth sample of signal sft) 
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The fundamental frequency (60 or 50 Ну) phasur magnitude and angle 
are then calculated as shown below: 
1 
MAGNITUDE = (REAL? + IMAG')* (peak value) 

ANGLE = arctan (TMAG/REAT) 


The DLP relays use both full-cycle and half-cycle "data windows" for 
the DFT calculations. "ће term "data window" refers to the time 
interval or, alternatively, the number of samples used in the ПЕТ 
summation. The actual DFT algorithm is recursive. А recursive ПЕТ 
uses a sliding data window where the calculations shown above are 
performed each time a new sample is available. The newest sample 
becomes k-N and what was previously sample К=1 is discarded. Such a 
recursive implementation maintains a constant phase angle reference 
ао that the phase angle does not rotate but remains fixed fur a 
steady sine wave input. Appendix II shows that the DFT calculations 
will produce the proper magnitude and phase angle when the waveform 
being sampled is a symmetrical sina wave. 


By definition, the DFT calculations shawn ahove extract the 
fundamental frequency component from the sampled waverorm. А 
full-cycle data window ПЕТ does а good but not perfect job of 
filtering out any dc offset ав indicated on page 1 of Appendix 111, 
but a haàlf-cycle data window DFT dees а poor job as indicated сп 
page 2 of Appendix 111. Appendix III uses the same current waveform 
from Appendix I. Equations (1) and (2) merely define the sampled 
values used by the ПЕТ and are not part ог іне ПЕТ calculations. It 
is the inability of the DFT to eliminate the dc offset that requires 
the use of the "software transactor" algorithm described previously. 


HARDWARE & ЕТЕМНАКГ - 


The DLP relays use three microprocessors; (2) 80C186 and (1) 
TM5320cC1n,. The HARDWARE DESCRIPTION and MODULES sections of each 
instruction book provides physical dimensions, front and rear view 
phoLographs, a module location diagram, an overall block diagram, 
and individual module block diagrams. 


The three modules, SSP, DAP, and БЕР, each contain a 
microprocessor and two EPROM sockets, When a firmware update is 
implemented the user will have to replace 1, 2, or 2 pairs of EPROMs 
depending upon the extent of the changes made to the software. 


Anti-Aliasing Filter - 


The ANI module contains the anti-aliasing filters. "Aliasing" 
іп an error related to sampling rate which is manifested by lower 
frequencies appearing in the sampled signal which are not present in 
the input signal.  Aliasing can be prevented by passing the input 
signal through an analog low-pass filter prior to sampling. This 
filter must have a cut-off frequency egual Lo or less than one-half 
of the sampling frequency, fa. For the DLP, Е ig equal to 16 X 
SYSFREG where SYSFREQ is either 50 Hz ог 60 Hz. Тһе lowest value of 
Га/2 is 400 Hz at 50 Hz, and the anti-aliasing low-pass filter ia n 
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three-pole butterworth design with a gain of -25 db at 400 Hz. This 
one filter design is used at both 50 and 50 Hz. 


Analog to Digital Converter - 


The DLP relays use one 12 bit analog-to-digital converter (ADC) 
with a multiplexed analog input system as shown in Figure 1. Since 
the total scan period for all the ADC inputs is approximately 65 
microseconds, the worst error due to sequential sampling is 1.4“ and 
"sample and held" circuits are not required. 


Six analog signals are sampled; (2) voltages proportional te 
VAG, VDG, and VCG and (3) currents proportional to IA, IB, and IC. 
The secondary voltages from the system PTs or CVTs are applied to 
auxiliary PTs in the DLP. These auxiliary PTs have а 14:1 ratio. 
The voltage samples are further attenuated 2:1 before they are 
applied to the Арс. Thea ADC's maximum input voltage is +5 volts dc, 
and each ADC bit represents; 


salt = 2.442 millivolts 
Referred to secondary system values ([(i.e., voltage applied to the 
relay), one bit represents 2.442*2*14 - 68.4 millivolts. The 
quantization error, q, is: 
а = 5*2 12 = 1.22 millivolts 


The quantization error referred to secondary system values, Qa, 15: 
dg = 1.22*2*14 = 34.2 millivolts 


Thus ап output of hexadecimal 001 (one bit) represents a secondary 
input voltage (at the terminals of the relay) between 34.2 and 102.6 
millivolts (peak АС). 


The secondary currents from the system CTs are applied to 
auxiliary CTs in the DLP. These auxiliary CTs have a resistor on 
their secondary to produce a vollage, labelled VIA, VIB, ог VIC in 
Figure 1, that 12 proportional to current. To achievo tha large 
dynamic range required for the current inputs while maintaining 
adequate resolution, the DLP uses two separate current inputs to the 
analog multiplexer. One input labelled "Ixió6" in Figure 1 covers 
the range of 0-20 amps peak and the other labelled "I" covers the 
range of 0-320 amps peak. Normally the "Ix16" input samples are 
uscd, but when an overflow on the “Tx16" input ів detected the 
time-equivalent sample from tha "I" input is used instead. "Ізі" 
simply means that the "I" samples are multiplied Ру 16. 


А secondary system current (і.е., current applicd to the relay) 
of either 20 ог 320 amps peak produces a 10 volt peak signal within 
the DLP. Like the voltage samples, the current samples are further 
attenuated 2:1 before they are applied to the ADC, 
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Figure 1 (Multiplexed Analog Input system] 


When "1:16" samples are used, as would be tha case for normal 
load flow, one ADC bit is: 


(5/211)424(20/10) = 9.77 milliamperes (peak АС) 


When "I" samples are used, as would be the case for a severe fault, 
one ADC bit 15: 


(5/211) *2*(320/10) = 156.25 milliamperes (peak АС) 


Considering the quantization error, one bit represente a secondary 
input current (at the terminals сЕ the relay} within the randes 
shown below: 


"Ix16" = [4,88 - 14.65) milliamperes {peak АС) 
"n - {78.13 = 234.38) milliamperes (реак АС) 


If the current applied to the relay should exceed 320 amperes peak, 
the resulting overflow is detected апа the current value is clamped 
at 320 amperes, 


Considering software calibration factors for the magnetics and 
АСС circuitry as well as software scaling factors, one bit on a 
voltage channel actually corresponds to 0.1 volts peak and one hit 
on an "Ix16" current channel actually corresponda to 0.01 amps peak. 


FUNCTIONS AND FEATURES 


Refer to the PRODUCT DESCRIPTION section of the instruction hook 
for the functions and features contained in a particular DLP relay. 


PHASE SELECTOR - 


Versions of the DLP are available for single-phase and 
three-phase tripping. For single-phase tripping, the phase-selector 
algorithm is used to determine which phase to trip. At the heart nf 
the phase-selector algorithm is fault type determination which is 
used for both single-phase and three-phase tripping to determine the 
fault type for the fault location algorithm. The fault component of 
positive-sequence current (Тір), negative-sequence current (Ia), and 
zero-sequence current (Та) are used to determine the fault type 
through a combination of level checks on Ia and Іп and the 
measurement of the phase angle between Іс and Тір. The fault 
component of positive-sequence current is ca culated ‘by Bubtracting 
the pre-fault current phasor from the fault current phasor. 


A test level for Ту, T2TVL, and a test level for Ig, IOLVL, за 
fixed in the DLP!s firmware. If |Ig| >= 1011, then a ground fault 
is declared and the angle by which Іс leads Гір is calculated to 
determine which of the six ground fault types is present. 


The theoretical phase angle criteria аге liated below: 


Ground Fault Angle (I, leads Тр) 

AG D" 

BG 120° 

CG 240" 

АВС ЕП" 

BCG 190" 

САС 300597 
"he actual criteria allows for а +/= error about the theoretical 
value. The phasor relationships between I, and Ij, for all 
unbalanced faults are shown graphically іп Figure 2. IÍ | Ipl < 


IOLVL, then a phase fault is declared and the angle by which 15 
leads Тір 15 calculated to determine which of three phase-to-phase 
faults might be present. If none of these angle checks (same values 
as for phase-to-phase ground faults} are satisfied, then a 
three-phase fault is declared provided that |Т.| < I2LVL. As a 
supplemental check, the fault type determined as described above 15 
compared agalnst "zone flags" which are generated by operation of 
the phase and ground distance functions. After the fault type iz 
determined, the phase-selector function produces four outputs, phase 
А, phage B, phase с, and multi-phase, which cause the DLP to trip 
the proper phase for single-line-to-ground faults and three phase 
for all multi-phase faults. 


AG Ва се 





ТЕ 
12 
АВ BC СА 
Ава BCG САО 
1 а 
2 180 
о UM SU 7 
во le МЕ заг г М 
А 


Figure 2 (Angle between Тук and I, for Unbalanced Faults} 
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FAULT LOCATION - 


Fault type determination is also required for fault location. 
Ггісг to running the fault location algorithm, the fault type must 
be known to select the proper input quantities. If fault type 
cannot be determined, fault location is not calculated, but the 
relay algorithms will still protect correctly. The fault location 
algorithm wlll be described and then the required input quantities 
based on fauit type will be listed. 


A "lumped-parameter differential equation" algorithm is used to 
implement a single-ended fault location estimate. "single-ended" 
means that only currents and voltages available at one end of the 
transmission line are used. Figure 3 shows a two terminal 
transmission linc with equivalent sources at each end. A 3-phaze 
fault exists at n*100 percent of the Linea length from terminal А. 
Assuming the DLP is located at terminal A, the voltage at А 15: 


VA = Леа ТА t Rp*(Ipa t Ів! 


where: n ~ fraction of line length from А to the fault 
51, = line impedance 
ТА = current at А 
Re = arc resistance at the fault 


ТЕРА = fault current at terminal А 
Тен = fault current at terminal Е 








Figure 3 (о Terminal Transmission Line) 
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Апу pre-fault load current is removed from the total current to 
produce Ігл and Igg. The term 2,,*T, can be rewritten as: 


dIa 
EL Ia = ЕТІ, + ГА at = Ши 





where: зт = R + jw*T. ғы = 24r*f} 
The term IZ i& a phasor corresponding to the sampled values iz that 
are created to implement the "software transactor" algorithm as 
previously discussed. The 12 samples are inputs to the ПЕТ and тя 
is the phasor output. 
The voltagc at À can be rewritten ag: 
МА = 17 + Ве? (Тед + L rp) 
A new term is defined as: 
Rra*Ifa = Re* (Тед + Тер) 
Irg 





where: Rfa = Rp*( 1 + 


— 


IfA 


to permit writing an equation for Уа in terms of quantities 
available at terminal А. 


МА = п"ја + ВЕКА ТЕРА 


The term Врд is а real number provided that Тед and ігр are in 
phase. For this to be true, the assumed system in Figure 3 must be 
homoqeneous. This means that the source and line impedances must 
have the same angle. While this is rarely the case in an actual 
power system, the assumption is made here to show that the fault 
location algorithm eliminates the effect of fault resistance if Efa 
is real. In practice, а typical non-homogeneocus system produces 
only a small error in the fault location estimate. 


Ма is now split into its real and imaginary components: 


VaA(real) Vs[ine] Va(imaginary) = VC[osins] 


VS + jVC = n*(IZS + IZC) + Rga^*(IgajS + ТЕЛО) 


үз 
ус 


n*IZS + КгА”ІғдӘ 
nh*IZzC + БА" ТАЈ 


Solving for "n" yields: 
VS*IpíC - VC*IgAS 


І25%ІғлАС - IZC*I pa 


proving that this approach eliminates the effect of fault resistance 
provided that the defined term ВЕд 1: assumed to be real. 


Like the apparent impedance seen by a distance relay, this fault 
location estimate 12 affected by current infeed from а tap or third 
terminal and Ру zero-sequence mutual impedance between parallel 
lines. Compensation for mutual effects using a portion of the 
zero-sequence current from the parallel line is not implemented in 
the DLP for either the fault location algorithm nor the distance 
functions. 


AS stated previously, the proper voltage, current, and IZ 
phasors musk be used in the equation for n, and these phasors аге 
determined by the fault type and faulted рһаве(в). For single-phase 
to ground faults, these phasors ага: 


a-g: Ұш МА 

ТЕ = (ТА) Е 

IZ = (ТА - Ig)*2,/POSANG + То“ Ко“ 21 /?ЕВАНМС 
b-g: У = Vp 

Ір = (ТВЈЕ 

IZ = (Ip - Ig)*Zi/POSANG + Ij*Kg9*Zi/ZERANG 
с=ч: V = Ve 

Је = (1с) р 

IZ = (Іс - Ig)*Zi/POSANG + І0%К9%21/2ЕБАНС 


For multi-phase faults these phasors are: 


a-b џ = VA - Ма 
ig = (1, - igip 
IZ = (Та = Ig) *2,/POSANG 
Те ^ (Тв - Ic)r 
12 = (Ір = Те) кај /РОЗАНО 
c-a V = Vo - Уа 
Te = (Те “ ТА)Е 
ТА = (Ic - Та) *2;/POSANG 
where: ( ік indicates that the pre-fault load flow current 
has been subtracted from the phase currents included 
in ( ) - 


OVERCURRENT/DIRECTIOWNAL FUNCTIONS - 


The overcurrent and diracticnal functions contained іп the СІР 
relays are listed below: 


FD - fault detector 


NT - negatbive-seqguence directional trip function 


NB - negative-sequence directional block function 
IPT - ground trip overcurrent 

IPB - ground block overcurrent 

IDT ~ ground overcurrent direct trin 

TOC - ground time overcurrent direct trip 

IT - trip supervision overcurrent 

тв - block supervision overcurrent 

PH4 — non-directional phase overcurrent direct trip 
Line Overload - (2) setting levels 

Та = line pickup overcurrent 


Fault Detector - 
The fault detector, FD, is a non-directional low-set overcurrent 


function with a factory set pickup level. Figure 4 is a block 
diagram for the FD function. 
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Figure 4 (Fault Detector Block Diagram) 
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The change (/\) in I I5, Of Ig magnitude is used together with 
Lhe total magnitude of Y ог Ір to determine FD operation. The 
AIl, |/АІЗІ, or ТОР quantity must exceed 0.2 amperes RMS or 
1-| or |та| must exceed 0.6 amperes RMS to produce an output. The 
1 quantity magnitudes are calculated each sample Бү comparing the 
respective IZ phasor output from the DFT with a memory value of IZ 
stered four cycles previously. Recall that all phase-current 
samples are converted to if samples to remove the dc offset and then 
fed to the ПЕТ algorithm which produces IZ phasors. Тћеве 
phase-current IZ phasors are then converted ta 18], IZ4, and Tig 
symmetrical component phasors. 


As an example, the positivessequence T74 magnitude calculated at 
the latest sample іє compared to the magnitude of IZ, g&tored in 
memory four cycles previously, and if tha absolute difference 


exceeds 0.2 amperes RMS, FD will produce an output. The four cycle 
"look back" comparison used for the /X quantities results in a four 
cycle "reset" time for the FD function. Since it is not only 


necessary that FD detect any system disturbance but also maintain а 
constant output during the disturbance, several inLernal DLP signals 
or flags are used to seal-in the FD output as shown in Figure 4. 
This level of detail із not shown on the functional logic diagrams 
in the instruction hooks. 


Ground Directional Functions - 


Figure 5 is а block diagram for the ground overcurrent and 
directional functions. Two negative-sequence dlrectional functions, 
NT ftrip} and НВ {block}, are used. It is well documented in the 
relay literature that negative=sequence directional functions are 
BEuperior to zero-s&Eenqguence current and/or voltage polarized 
directional functions particularly when zero-seguence mutual 
coupling is present between parallel lines. ні supervises the ірі 
pilot-trip function and provides selectable directional control 
(1.е., torque control) of the direct-trip backup functions, IDT (50) 
and TOC (51), HB supervises the IPB pilot-block function. 


NT and NB are implemented using an "amplitude comparator" 
measurement. Ап amplitude comparator is а generic type of relay 
function measurement technique that compares the magnitude of an 
OPERATE quantity versus the magnitude of a RESTRAINT quantity to 


determine if the relay function should operate or not. If the 
OPERATE quantity is larger than the RESTRAINT quantity, thon the 
relay function operates, For НТ, the simplest amplitude comparator 


design is: 
OPERATE = [у = I5%Zp| 


RESTRAINT = |V, + Т-%7р| 


where: In = negative-seqguence current at relay 
V4, = negative-sequence voltage at relay 
2р = relay "reach" impedance = |Zp| /POSANG 


=13- 
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Figure 5 (Ground Directional and Gvercurrent Functions) 
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However, for some system conditions, the negative-sequence voltage 
at the relay may approach zero. With V4 approaching zero, NT would 
be near it's balance point, and would not operate reliably for an 
unbalanced fault in the forward direction. Consequently, a 
compensation technique is used which creates а reliahle V; 
polarizing quantity in the relay even though Vj, at the relay may be 
zero. This compensation modifies the OPERATE and RESTRAINT 
quantities as shown below: 


OPERATE = |V; = (1+) +1.%2р| 
RESTRAINT = |У; + (1-K)*15*Zp| 
where: К = 0.05 = offset compensation 


For NB, the OPERATE and RESTRAINT quantities are simply 
interchanged. 


The actual DLP OPERATE and RESTRATNT quantities are tabulated 
helow; 


НТ NB 
OPERATE [v5 - 1.05*I5*Zpl |2*V4 + 2*T5*Zp| 
RESTRAINT |V5 + токл| |2.05*®1„*2р| 


where: Ар 
ZR 


20 /POSANG for а 5 ampere DLP 
100 /POSANG for а 1 ampere DLP 


For NT to operate, OPERATE must exceed RESTRAINT by 1.0 secondary 
volta RMS. For HB to operate, OPERATE must exceed RESTRAINT by 0.5 
secondary vults RMS. The NT RESTRAINT guantity is КЕ + Трка 
rather than the theoretical | Мо - (1-К)%1-%7;|. The elimination of 
the (1-K) factor creates a larger magnitude restraint quantity for 
increased security. For the КВ function, both the OPERATE and 
RESTRAINT quantities are essentially doubled in magnitude compared 
to the corresponding НТ values. А blocking function at one end of а 
protected line should operate ав fast as or faster than the tripping 
function at the other end for a through fault condition. Since the 
DLP employs an enhanced amplitude comparator called an "energy 
comparator", doubling the NB quantities compared to NT assures more 
enarqy in NH and consequently a faster operate time compared to NT 
at the other end. 


The basic advantage of the energy comparator іє that an 
integrating algorithm replaces the magnitude comparison of the 
simple amplitude comparator. ІТ the difference between the OPERATE 
and RESTRAINT signal 1з positive, then the integrating algorithm 


integrates "up" toward the operate threshold level. If this 
difference is negative, then the integrating algorithm integrates 
"down" away from the operate threshold level. This allows for 
modifying the response of NT and NB by adding restraint "energy" or 
bias when the fault detector, FD, is not operated. This is 


indicated in Figure 5 by the BIAS ARTIST inputs to NT and HB. When 


FD operates indicating a system fault or disturbance, this bias 
energy is removed, and should NT or NB begin to integrate "up" 
toward operation the integration must overcome the initial restraint 
energy. Simply stated, at the instance a fault occurs both WT and 
МВ are biased toward non-operation to increase security. 


Other Lhan the factor of two mentioned above there are other 
differences in the NB OPERATE and RESTRAINT quantities compared to 
the theoretical quantities. The RESTRAINT quantity is |2.08#1) #2 | 
rather than |2*V4 - 2%(1%к)%1-%2р|. Elimination of 2*V, creafes а 
smaller RESTRAINT quantity which increases dependability and speed. 
The OPERATE quantity із |2*V4, : 2*1I5,*Zg| rather than |2%У- + 
2%(1-К)%І-%7-|. Elimination of the (1-К) factor creates a larger 
OPERATE quantity which increases dependability and speed. 


Ground Pilot Overcurrent Functions - 


IPT ANDed with NT is the ground directional-overcurrent pilot 
trip Yunction, and IPB ANDed with WA is the ground directional 
overcurrent pilot block function. These functions aro labelled TRIT 
and BLOCK in Figure 5, and they are only in service if SELZ2U = D or 
^ апа one nf the pilot schemes (i.e., POTT, PUTT, BLOCK, or HYBRID) 
has been selected. Тһе IPT and IPB operating quantities аге; 


IPT: 3*|rg] - зекта|т. | 
IPB! 3*|Ig| - 3*kBx|1, | 
неге: KT = 0 nr 0.1 
ЕВ ~ 0.066 

For the BLOCK and HYBRID schemes, KT = 0.1. Fur Lhe POTT and PUTT 
Schemes, KT = 0. The КТ and KB constants are not directly user 
selecLable but are set to the appropriate value based on the SELSCHM 
setting. Figure 5 shows that the local NT, zone 1 distance, and 


zone 2 distance functions (via OR 103) block NB*IPB at AND 2. This 
local "trip over block" preference is part of the transient blocking 
logic and is not shown on the functional logic diagrams in the 
instruction books, 


IPT and ТРА utilize positive-sequence current restraint, 
J*KB*HI,], to Limit their reach for external faults and to prevent 
possible operation due to system or load dissymmetry under maximum 
load flow conditions. This feature increases scheme security when а 
BLOCKING or IIYBEID scheme is selected. When a POTT ог BUTT tripping 
scheme is selected, the positive-sequence current regtraint is 
eliminated from the NT function (ЕТ = 0) to assure that it operates 
for all remote terminal faults under the most severe system and load 
conditions. 


Ground Backup Overcurrent Functions - 
An instantaneous (ТОТ) and time overcurrent (TOC) function 


provide direct-trip backup protection for ground faults. The ТОТ 
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operating quantity is: 
з*|1] - 3*ED*iIjl 
where: HD = 0 or 0.34 


In some DLP relays, KD is fixed at 0.3. Tn other DLP relays, KD can 
ре set to 0 or 0.3. The setting determination is more difficult 
(see the instruction books) when KD - 0.3, but overreacn 
possibilities are greatly reduced while sensitivity is increased 
after the remote breaker trips. 


The TOC operating quantity is simply 3*|Ij|. Depending on the 
version of the DLP relay, the characteristic curve shape is either 
fixed or user selectabla. Both ТОТ апа ТОС can be separately 
controlled by the NT directional function or operated as 
non-directional functions. This "directional control" is analogous 
to "torque control" in EM relays since the ТОТ or TOC algorithm does 
not begin executing until КТ operates. 


FPhase-Current Actuated Overcurrent Functions - 


IT (trip) and IB (block) are two low-set overcurrent supervision 
functions each of which have seven outputs. às shown in Figure 6, 
three outputs ЕТТА, ITS, ITC} are per phate, three outpute (ІТАВ, 
ITBC, ITCA) are per phase-pair, and one (IT) is a logical OR of the 
per phase outputs. ‘lhe settings, FUIT and ЕСІН, are in terms of EMS 
secondary phase current, As shown in Fiqure 6, the threshold level 
for the per phase-pair outputs is /43*FUIT. 

IB provides overcurrent supervision of the zono 4 distance 
functions. IT provides overcurrent supervision of all other 
distance functions {refer to the next section). Ground distance 
functions are supervised by ITA, ITB, or ITC. Phase distance 
functions are supervised by ITAB, ІТЕС, or ІТСА. IT and TB are also 
used for other purposes as shown on the functional logie diagrams in 
the instruction books. Supervision of the distance functions by 
these two overcurrent functions is primarily intended to provide a 
fast reset of the distance functions; it does pot provide blown fuse 
protection since PUIT and PUIB are normally set well below the 
full-load current. À separate function is provided in the DLP 
relays to detect a blown AC potential fuse. 


ЕНА із an instantaneous non-direcbional overcurrent direct trip 
function. It is intended to provide direct tripping for multi-phase 
faults, and it operates on the highest of the three delta currents, 
la 7 Ір, Ip 7 le, or Те ~ I}. The pickup setting, PUPH4, іє based 
on the delta current magnitude not the phase current magnitude. 


There are two overcurrent functions (two settings - PULVI and 


PULV2) associated with the Line Overload function. Each function 
operates on the highest of the three phase currents, Та, Тв, or Іс. 
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Figure в (Overcurrent Supervision Functions) 
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Line Pickup Overcurrent Function - 


The overcurrant function, Il, associated with Line Pickup 
operates on the magnitude of positive-sequence current. Operation 
of Line Pickup is described in the instruction books. 


DISTANCE FUNCTIONS - 


1n this section, the phasor inputs to the Phase Angle Comparatcr 
for each distance funetion are listed. This information will allow 
the user to analytically determine the steady-state response of any 
distance function for any fault. Ап example of such an analysis ік 
given at the end of this section. 


A Phase Angle Comparator is a measurement technique that 
determines the amount of coincidence between two or more phasor 
inputs (іп electrical degrees) to determine if the relay function 
should operate or not. References 1, 2, and 3 provide more detailed 
information about Phase Angle Comparators. In the DLF, the 
associated phasor inputs (magnitude and angle) are first calculated 
and stored in memory then their electrical coincidence is determined 
using a simple software routine. 


For each zone of distance protection there are six (6) separate 
distance functions. There are three phase distanca functions - one 
per phase pair (1,е., AB, BC, and СА), There are three ground 
distance functions - one per phase [1.e., А, B, and с). 


The following definitions pertain to all of the distance 
functions: 


In = phase А current at relay 
Ig = phase B current at relay 
Та = phase C current at relay 


Ig — zero-sequence current at relay 

Ма — phase А to ground voltage at relay 
Vg = phase B to ground voltage at relay 
Ve — phase C to ground voltage at relay 


( )1 = positive-sequence component of | ) 
( ја ^ neqative-sequence component of | ) 
( ум = memory fpre-fault) value of | ) 


ZXl ^ zone x ров.-вей, reach setting - ZxR /POSANG 
2x0 ~ zone x zero-seg. reach setting - ZxR /ZERANG 


[note: the magnitudes of Zx1 and 2х0 are identical - these 
quantities differ in phase angle only] 


240R = zone 4 offset reach multiplier 
КО1 = zone 1 zero-seguence compensation factor = “ак 
КО = gero-sequence compensation factor 


In the DLP relays, КО is equal to гр of the line divided by 21 af 
the line. 
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Fone 1 Distance Functions = 


The variable МНО zone 1 distance Phase Angle Comparators 
fixed at 90° coincidence to produce a circular characteristic. 


Phase Distance - variable MHO: 


АБ: 


BC; 


СА: 


(ТА - Ір) *211 - (Và - Ур) 
(Ма 7 Vadim 

{ТА — Ів) “211 

(Ip - Іс) *211 - (Vg “ Ус) 
(Ve 7 Ус) 1м 

(Ig - ІС) #211 


(Ус = Уд)1м 
(Ie = In} *211 


Ground Distance - variable МПС: 


А: 


(I, - То) #211 + К01%10%210 
(Vadim 

ТА) 5*Z11 

19%211 


{Тв - Іу)%211 + KOl*Ij*Z10 


(Vn) 1M 
(Ір)2%211 


Ip*211 


(Іс - То! 5211 + ЕО1#Т 7210 
(Ус) 1м 

(Ic)2*211 

то #211 


Ground Distance - reactance: 


А: 


(ТА - 10) *211 + KO1*Ty*Z10 
(Та)2%411 
1211 


(Ig - Ig)*Zll + К01%1,%210 
(Ің)з%211 
Tye al 1 


(Іс - Ig} *Z11 + К01%10%210 
(Те) #211 
104211 


Ус 


are 


The supervising MHO that is part af the total reactance function 
has the same phasor inputs as the zone 1 ground MHO unit except that 
211 is determined by the setting 71617. 


load flow conditions. 


in the instruction books. 


The 2150 reach adapts to 


Reference the CALCULATION OF SETTINGS section 


zone 2 Distance Functions - 


The number of electrical degrees of coincidence required to 
produce an output from the Phase Angle Comparators of the variable 
МНО zone 2 distance functions is determined by the settings Z2PANG 


and #2 САН с. 


Phase Distance = variable МНО: 


АВ: 


BC: 


CA: 


(ТА - Ig)*Z221 - (Уд - Уа) 
(МА 7 Ув 1M 

(ТА - Ір)%221 

(Ig - Іс)%221 - (Vg - Vc) 
(Vg — Ус) 1м 

(Ig - 1с) #221 

(Tp = TA) *Z21 = (Ул = Va) 
(Vc = Vadim 

(Іс - Тд)*721 


Ground Distance - variable MHo: 


А: 


(Та - 16) +221 + ЕК041,%220 
(УАЗ ам 

(ТА) 3221 

Та? 221 


(Ін ~ Ig)*Z21 + ко» “220 
(Vn) 1M 

(Ір) >*Z21 

ЕТЕ 


(Ip - Ig)*Z21 + E0*IQ*Z20 
(Voi 1M 

(10)2%721 

125721 


"one 3 Distance Functions - 


The number of clectrical degrees 
produce an output from the Phase Angle Comparators of Lhe variable 
МНО zone 3 distance functions is determined by the settings Z3PANC 


and Z3GANG. 


Phase Distance - variable МНО: 


AB: 


BC: 


(Та - Ig)*231 - (Và - Ур) 
(Ма 7 Ув)ум 
Пат ig) *Z21 


(Ip - Іс)%231 - (Vp - Ус) 


(Vg – Уе) ум 
(Ip - Іс) +231 


of coincidence required to 


СА: (la - 1д}*731 - (Ме - Уа) 
ivo > VAIM 
Те - ТА *911 


дъулита Drsbainos - variable MEC: 


а: iI, - То) *231 + KO*Ig*Z30 ТА 
МАЈ 1м 
I UPS 
Tk hal 


м 


1 


5: lip - 1117843. + Kuüu*I4*zZ3u “ы 
(Ур!) 
Тр! 5431 
Ty tha 
C i г - Ins #2 31 + Кож РКЕ - Vo 
M 
E 


^ 
1 


= н 


*(]:172 


E 
7a 
1 


7 


Ld 


| 1 М 
d 2 1 
14721 


га 


Fone 4 Готен Functions - 


Тие qurber of eliectricai deqrecs ot coincidence required то 
стосе an cctput frem the Thane Angle Comparators of the variable 
MEG zere 4 distance functione if determined by the settings Z4PANG 


ага УЗАЛА. 


Fiese Distarce - variable МПО: 

AD: ila Ig; #441 - (МА - Мы) 
“ТА - Газа ХАЧОКААЧЈ - (va - Voli 
Ід - Тај 241 

BU: fap T сојева (Мы Yol 
д rol: tAÉAORY*ZA]l - Ма - Уң!-м 
ГА - Те кад 

CA: ile- пума - We- Ya 
[la - Та) 7 *260R* ЛАТ - iVa - ҮБ 1M 


те ТДд}*Й41 


vhe abovo phasor inputs are for the forward reach cage where 
ата; = 0. For the reverse reach case, where SELZ4D ~ 1, simply 
иш à minas sign in front of each питеспе exprEesslor, 


iroand Gistarce - varijabie МНО: 
A: "ТА Та 3291 + Кот "а - Ул, 
“МАЛ ТМ 
Inia Adal 
Tatra 


-22- 


B: (Ір - Ig)*Z41 + КО%1,%740 - Vp 
(УВ) 1M 
(Тр) „Ха41 
1 2241 


C: {те - Ig)*Z41 + коктока40 - Vo 
(Vc) iM 
(16) 2^241 
1ogf*£41 


The above phasor inputs are for the forward reach case where 
SELZ4D = 0. For the reverse reach case, where SELZ4D = 1, simply 
put a minus sign in front of each current expression. 


Analysis Using Phasor Input Coincidence - 


The following is an example of how calculation of phasor input 
coincidence can be used Ко determine if а distance function should 
or ehould not operate for a given fault condition. Гог tha sample 
transmission system of Figure CS-1 in the instruction books, assume 
а DLF relay is located at terminal Able with а bolted A-G fault 
located at terminal Delta. The phase А, воле 3, ground MHO function 
will he analyzed. The instruction book calculation shows that the 
apparent impedance seen at terminal Able for this fault is 15 /80.31" 
secondary ohms assuming no pre-fault load flow (1.е., voltages El, 
Е2, and ЕЗ are all at 0°). For this condition, the phase A, zone 3, 
ground МНО function is at a virtual balance point with Z3CR (221) = 


16.00. The evaluation below takes into account pre-fault load [Clow 
(1.с., ДЕТ - O°, /E2 - -35', ЈЕЗ - -30"). 

ЗБЕ (231) = 16.00 

POSANG = 65° 

ZERANG = ЈА" 

KO = 3.2 

Z3GANG ~ 90 


The fault currents and voltages at the relay are listed below: 


pre-fault Уд = 64.87 „==> 6" 


Va = 59.93 /-3.9* 
(Va) = 62.82 /-4,9' 
ТА = 3.915 /-53.9 


(Та) = 1.099 /-111.2: 
Ір = 0.493 /-105.1" 


The four phasor inputs are calculated to be: 
1) (ТА - Ij)*Z31 + КОЗТа? 210 ~ Ул = 27.26 /72,8" 
2) (Valin = 64.87 /-5,6* 


3) (15) 1231 = 17.58 /-26.2* 
4) 10*231 = 7.89 /-20,1' 


The number of electrical degrees of coincidence, considering all 
four phasors simultaneously, must be determined. Іп this case, all 
four phaser inputs must overlap for at least 90%. This overlap ог 
caincidence is equal to: 


180" + [-26.2") - 72.8" = 81" 


Since 90% af coincidence is required for operation (42GANG = 90), 
the phase A, zone 3, ground МНО function will not operate for this 
fault condition. 


This method of determining distance function operation is 
preferable to plotting impedances on an H-X diagram since the 
variable MHO characteristic changes its size and position on the R-X 
diagram depending on pre-fault lead flow, system impedances, and 
fault type. Іп this example, a "no pre-fault load flow" calculation 
indicates that a 16.00 ohm reach setting is just barely adequate. 
However, when pre-fault lead flow is taken into account the 16.00 
ohm reach is not adequate. Іп fact, the reach must be increased to 
approximately 18.00 ohms to assure operation for this condition. 


Current sensitivity - 


Any distance function requires a minimum value of current to 
assure operation. In the instruction Looks for the DLP relays, 
formulas for determining distance function current sensitivity arce 
located near the end of the CALCULATION OF SETTINGS tab section. 
These formulas are based on the inherent sensitivity of the distance 
functions and do not take inte account the pickup settings of the IT 
or IB overcurrent supervision functions. 


There is no one value of fault (or test) current helow which а 
distance function will not operate and above which it will reach 
100% cf its set reach. Instead, the actual reach decreases compared 
to the set reach as the current decreases. For electromechanical 
distance relays reach versus current was plotted as a "bullet 
curve." This decrease in reach is often expressed as "pull-back." 
Fur instance, if the set reach is 4 ohms but the actual renuch at а 
certain level of current is 3.6 ohms, then the reach "pull-back" is 
0.4 ohms or 10% of the set reach. Тһс current sensitivity formulas 
in the DLP instruction books utilize this concept of reach 
"ръ11-Баск. " 


Polarizing Voltage Sensitivity - 


The МНО distance functions utilize positive-sequence voltage 
with pre-fault memory for the polarizing quantity. Тһів means that 
for any bolted zero-voltage fault, other than a three-phase fault, 
there will Бе sufficient polarizing voltage to keep the distance 
functions operated as long as the fault persists (assuming 
gEutfficient fault current). For a bolted three-phase fault, Lhe 
phase distance functions will operate transiently because of the 
pre-fault memory voltage hut will eventually reset if the memory 
voltage resets beforc tho fault іс cleared. 
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The minimum value of positive-sequence polarizinq voltage 
required at the DLP to keep the МНО distance functions (all zones 
operated for as leng аз a three-phase fault persists іс 5% of 6? 
volts or 3.35 volts RMS. This fact is important when providing zone 
2 backup protection on а "short" line. Depending on the source to 
line ratio, the voltage at the relay may approach zero for a 
three-phase fault at the remote terminal. Unless the relay voltage 
is 3.35 volts or greater, the zone 2 function will not stay picked 
up leng enough to time out the zone 2 timer. 


It should be noted that some level of polarizing voltage 
sensitivity ір desirable. For a three-phase fault with arc 
resistance directly behind the relay, the effect of pre-fault load 
flow in the non-tripping direction causes the arc voltage to shift 
in phase towards the operating région of the distance funotion. The 
effect of polarizing voltage sensitivity is to raise the threshold 
of the polarizing voltage above that produced by the arc drop and so 
prevent operation. 


Operating Time Curves - 


сте of the phasor inpuls listed above that is commen to all of 
the distance functions is "Iz-V" which is often termed the operating 
quantity.  IZ-V has two basic forms that differ between phase and 
ground distance functions as shown below. 





IZ-V 
phase variable MHO: Ilop*?Ri 7 "ag 
ground variable МНО: (Із-І)%?2рі) + КОЯТО? ВО 7 Vy 
where: Тар = phase-to-phase current at relay (e.g., I, - Та! 
Ig — phase current at relay  (e.q., 14) 
Ту = zero-sequence current at relay 
Vag = phase-to-phase voltage at relay (в.п., Va - Уд) 
Уа = Phase-to-ground voltage at relay e.g., Уд) 
Zp, = positive-sequence reach setting - Znil /POSANG 
Әр = Zero-sequence reach setting = Тара / ZERANG 
КО = zero-sequence compensation factor 


The operating times of the MHO distance functions are directly 
related to the magnitude of the operating quantity. Figures 7 and 8 
plot operating time versus the magnitude of IZ-V. These two figuras 
ага only applicable to DLP revision C relays. The inverse nature of 
these time versus IZ-V curves results from an intentional design 
effort to optimize speed versus security. Simply stated, a heavy 
cloge-in fault (high 14-V) has the fastest tripping time berause 
such a fault poses the greatest threat to power system security. 
Power system security is impacted less for remote line-end faults 
(low IZ-Vj, but relay security is more critical because the fault is 
closer to the halance point of the distance function. 1t is 
thercfore prudent to increase the tripping time for low values of 
IZ-V to assure relay security without sacrificing power system 
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Figure 7 (Operating Times with PT Potential Sources) 
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Figure 8 (Operating Times with CVT Potential Sources) 
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security. Stated another way, puwer system security could be 
severely compromised if a zone 1 distance function were to 
over-reach for a remote bus fault. 


For the DLP revision с relays, the setting SELPRIM can bc set to 
either ü(CVT PRI), 1(СУТ SEC), 2 (РТ РЕТ), ог ЗЕРТ SEC]. This one 
setting determines two different aepects of the DLF's operation. 
Firet, the setting determines how PRESENT VALUES are displayed. 
With SELPRIM = О or 2, the PRESENT VALUES (currents, voltages, 


watts, and vars} are displayed and stored as primary valucs. With 
SELPRIM = 1 or 3, the PRESENT VATIES are displayed and stored as 
secondary valuez. All settings are expressed in terms of secondary 


values, regardless of the SELPRIM setting. 


Secondly, the setting determines the amount of filtering used in 
the DLP's distance functions ta overcome transient error signals 
that may Бе present in the AC voltages. When magnetic voltage 
LEransformers (PTs) are weed, SELPRIM = 2 or 3 should ре selected, 
depending upon whether PRESENT VALUES are to be displayed in primary 
or secondary values. When capacitive voltage transformers (СТ) 
are used, SELPRIM = 0 or 1 should be selected, depending upon 
whether PRESENT VALUES are to be displayed in primary ог secondary 
values. With SELPRIM = D of 1 the operating time of the distance 
functione will be &lower at lower valuee of operating signal (13-1) 
where the transient error signals associated with СУТя become 
significant. Note that when CVTS are used, a setting of SELPRIM = 2 
or 3 may result in zone 1 overreach for line-end faults. Therefore 
SELTRIM must Бе set to 0 ог 1 when CVTs are usi. 


Figure 7 applies when SELPRIM — 2 or 3, and Figure 8 applies 
when SELPRIM = U or 1. Note that IZ-V is expressed as a per unit 
value. The per unit "base" is 67 volts for a ground distance 
function and /3%57 = 116 volts for a phase distance function. To 
use these curves, thc appropriate currents and voltages at the relay 
must be determined from а short-circuit analysis and the DLP relay 
settings must be selected for the particular zone сі protection 
being analyzed. Figures 7 and в are valid for all the variable МПО 
distance functions regardless of the zone (i.e., zone 1, zaone 2, 
zone 3, or zone 4). 


The following іс an example of how to calculate operating time 
using Figures 7 and 8. For the sample transmission system of Figure 
Се-1 in the instruction books, the operating time of the zone 1 
ground and phase distance МПО functions will be computed for а DLF 
revision г relay located at terminal Able. The zone 1 reach ін 905 
of the line А-В positive-sequence impedance and the faults are 
located at 60% of the line А-В positive-seqguence impedance in front 
uf the relay. 


he relay settings are: 


ТЛЕ = 5,41 
215Е - 5.40 


2160 = 
POSANG 


SERANG 


B^ 
74 


Е И са 
e 


Therefore: 


5.40 /B5"* 
R.4f /74% 


Zn] 
ва 


For ап a-g fault: 


Ја = 8.6H6 /—73.4" 
Та = 2.696 /-95.7° 
МА = 50.06 /-1.4: 
І2-У — (I&4-Ip)*Zpg; + KO0*Ig*Zpg - Уд - 22.3 /-1.637 
22.3 
Iz-V [per unit) = — = 0.333 
67 


From Figure 7, the operating time is 21.5 to 24 milliseconds. 
For a B-C fault: 


Ің-Іс = 20.53 /-175" 
Ұқ-Ус = 73.94 у-9о0“ 


1Z-V = (Ig-Ic)*Zpi - (Vpg-Vc) = 36.97 7-807 
36.97 


IZ-V (рег unit} = = 0.319 
116 





From Figure 7, the operating time is 22 to 24.5 milliseconds. 


For this example, Fiqure 7 or Figure 8 will produce the same 
operating times since the curves overlap when IZ-V is greater than 
approximately 9.2 par unit. Note that the operating times from 
Figures 7 and 8 are based on SCR trip outputs. If contact outputs 
are used, 4 milliseconds must be added to the operating times. 
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APPENDIX I 


General equation for fault current (including dc offset): 


-tL 
i(t) = I*sin(ut + а - 6) - I*sin(a - G)*c 
B = arctan (WL/R) = impedance angle of power system е ;— 85 
ш := 2-77-60 (à = 376.991 
и = angle on voltage wave at which fault occurs g i= -5 
" | 
L/R = Т (Гог power system) T = Бал |85: — 
190 ш 
т 
Let: ф t= (m - 8): — ФЕ =1.571 radians (ф = -90 degrees) 
180 
t = 0,1 ..657 


note: w and 1/Т are multiplied by 0.0001 to allow time range above 
to be expresse in seconds 


à r= ш- „0001 T := T. 10000 
ш = D.038 T — 303,192 
“ 
Ti 
i(t) зе siní(u-t + 4) - sin($).e 


1 (+) 


0 -4 060? 
t-10 


(SECONDS) 
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APPENDIX I 


Primary Impedance - В + 701. Replica ішрейапсе = Rz + JOL? 


where: L/R - L2/R2 = T H2 t= 1 La := T 
( AA = 8 = ВА) 
d 
v(L) :— iit) R2 + Е) L2 
dt 


vit) — voltage across replica impedance 


i(t) зн i[t):109 ( ift) is multiplied by 10 to show up better 
on the plot below ) 


20 





wD pe, . . 
Ü -4 . 0667 
t.10 
(SECONDS } 
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APPENDIX II 


DISCRETE FOURIER TRANSFORM 


| (equation to define 
у := 1,523 2'п.— + — sampled values] 








2 2 IT] 180 
PEAK := JR + I ANGLE те atan —! -—— 
Ri т 


= al 


РЕАК - 1.5 ANGLE = 45 


APPENDIX III 


FULL-CYCLE DISCRETE FOURIER TRANSFORM 


H :— 16 k := 0,1 ..15 Т іш D.0303191 
| Kk 
k т | T N-60- T. 
(1) i = яігп|2-л- - 90- - sin,-8D. е 
К N 180 | 180 


( ароуе equation defines sampled vaiues for current with 
dc offset - see Appendix І } 


i ( 1 = sampled values of current ) 





g Е 15 


real (Е) and imaginary (I) equations for ПЕТ: 


à k 2 ка | k 
E := – > і .віп|2-Я.- | ќ= -- > і -сов|2-п.- 
ме c N Н — К қ. 
К 
Е = 0,132 І = -0.262 
2 2 I| 180 
PEAK :- Е + I ANGLE i= абап|- | · — 
в) т 
PEAK = 0,971 ANGLE - -82.194 
actual peak — 1.0 actual angle = -85 
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APPENDIX III 


HALF-CYCLE DISCRETE FOURIER TRANSFORM 





N i= 8 К: 9,1.,7 T :- 0.0302191 
k - 
4 л т N-2-.60-T 
(2) 1 = 5in|-.— - 50.—— - suin|-560:—— :е 
к N 180. 180. 


( above aquation defines sampled values for current with 
dc offset - see Appendix I ) 


( 1 = sampled values of current ) 





0 К 7 


2 (Kk 2 | k] 
E := —. 1 "аа m-— І := LN 1 -cosl m-m-—;, 
М —— kK (ON N — k м: 
К k 
Е = 1.097 T = -0.87 
2 2 I 180 
РЕАК :— ЈЕ + I ANGLE ‘= atan|- ---- 
РЕАК — 1.282 ANGLE = -31.343 
acLual peak = 1.0 acLual angle = -85 
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